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I. INTRODUCTION 



t-H , It is generally believed that only massive fermions have Bohmian trajectories but bosons do not. This is usually 
attributed to the impossibility of constructing a relativistic quantum mechanics of bosons with a conserved four- 
vector probability current density with a positive definite time component. However, it has now been shown [jlj that 
a consistent relativistic quantum mechanics of spin and spin 1 bosons can be developed using the Kemmer equation 
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(inp^d* + m c)ip = (1) 

where the matrices (3 satisfy the algebra 

P„ Pv P\ + P\ Pv (3fi = P„ g vX + (3 X . (2) 



, The 5x5 dimensional representation of these matrices describes spin bosons and the 10 x 10 dimensional represen- 
tation describes spin 1 bosons. The fact that a conserved four-vector current with a positive definite time component 
OA , can be defined using this formalism can be seen as follows. Multiplying ([!]) by Pq, one obtains the Schrodinger form 
■ of the equation 



dip 

ih— = [-ihcfiidi - m c 2 /3 ]ip (3) 



dt 

where pi = Pq Pi — Pi Pq. Multiplying (|l|) by 1 — /3 2 , one obtains the first class constraint 



Cd z?i/3 I /3 2 9^ = -moc(l-/3 2 )^. (4) 

It implies the four conditions A = V x B and \7.E = in the spin-1 case. The reader is referred to Ref. (j]] for further 
^ ■ discussions regarding the significance of this constraint. 

If one multiplies equation (|J) by if>^ from the left, its hermitian conjugate by ijj from the right and adds the resultant 
■ equations, one obtains the continuity equation 

a ' 



dt 

This can be written in the form 



+ di^pii> = 0. (5) 



0" = (6) 



where — — G^a" (with a v a u = 1 where a v is the unit four- velocity of the observer), = —m,Q(P'tl){P^P v + P v p^ — 
g^ip is the symmetric energy-momentum tensor so that Goo = — mo c 2 ip' if) < 0. Notice that s M s M = G^G^ > 0, 
so that is time-like. Thus, it is possible to define a wave function <fi — \/ mo c 2 / E ip (with E = — j Goo dV ) such 
that (jr (f> is non-negative and normalized and can be interpreted as a probability density. The conserved probability 
current density is s M = —Q^o/E = (0t <f>, Pi <j>) ffl. 

Notice that according to the equation of motion (||), the velocity operator for massive bosons is cfii, so that the 
Bohmian 3- velocity can be defined by 

dXi _ x Ui Si '0 t Pi^ rn\ 

v i = ~~TT = 7 u l =c— = c— = c — — — . (7) 

dt uq so tpi ip 
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It follows from equation (3) that c/% is the velocity operator whose eigenvalues are ±c. Therefore, v^v^ — 0, and 
so the Bohmian velocity is always timelike. Integrating equation (7), one obtains a system of Bohmian trajectories 
xi (t) corresponding to different initial positions of the particle. In Bohmian mechanics one assumes that the initial 
distribution of the positions is given by ^(O)! 2 . The continuity equation (^) then guarantees that the distribution 
will agree with quantum mechanics at all future times. The (Gibbs) ensemble averages of all dynamical variables in 
Bohmian mechanics will therefore always agree with the expectation values of the corresponding hemitian operators 
in quantum mechanics. 

The theory of massless spin and spin 1 bosons cannot be obtained simply by taking the limit mo going to zero. 
One has to start with the equation S 

ihp^^jj + m a cT^ = (8) 
where T is a matrix that satisfies the following conditions: 

r 2 = T (9) 

r^ + ^r = /v (io) 

Multiplying (||) from the left by 1 — T, one obtains 

/3„0"(I» = O. (11) 
Multiplying (||) from the left by d\ (3 X j3 v , one also obtains 

d x t3 x (3 v (TiP) =d v {Til>). (12) 

It follows from Q and @ that 

□ (Tip) = (13) 
which shows that T ip describes massless bosons. The Schrodinger form of the equation 

=-iKchdi(T^) (14) 

and the associated first class constraint 

ihPiffi d^ + m c(l- (3l)T^ = (15) 

follow by multiplying (||) by (3q and 1 — 0$ respectively. Equation (14) implies the Maxwell equations curli? = 
-(n/c)d t H and curLff = (e/c)d t E if 



T^ T = (1/ Vm^)(-D x ,-D y ,~D z ,B x , B y , B z , 0,0,0,0) (16) 
The constraint (15) implies the relations divE — and B = curL4. The symmetrical energy-momentum tensor is 



and so the energy density 
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e„» = -^<KAA + - 9^)^ (i7) 



£ = -e 00 - T ^ = \ &e + B-B] (18) 



is positive definite. The rest of the arguments are analogous to the massive case. 
The Bohmian 3-velocity for massless bosons can be defined by 

Using arguments similar to the case of massive bosons, it is easy to see that the Bohmian velocities for massless bosons 
are also timelike. Integrating equation (19) with different initial positions, one gets a system of Bohmian trajectories 
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for the photon. Neutral massless vector bosons are very special in quantum mechanics. Their wave function is real, 
and so their charge current = <fi T j3^ <j> vanishes. However, their probability current density does not vanish. 
Furthermore, s, turns out to be proportional to the Poynting vector, as it should. 

In this paper we compute Bohmian trajectories for photons for certain simple but interesting cases. Integral 
curves of the Poynting vector for localized wave packets in classical electrodynamics were first plotted by Prosser [Q . 
They are lines of energy flow in classical electrodynamics and cannot be interpreted as particle trajectories. A particle 
trajectory interpretation of these curves is possible only within the context of a proper relativistic quantum mechanics 
of indivisible photons. This is what we have done to calculate Bohmian velocities and hence trajectories for photons, 
carrying the entire interpretational package of Bohmian mechanics. It is only incidental that such trajectories happen 
to coincide with the integral curves of the Poynting vector for single photons. However, in the case of two photons, the 
Bohmian trajectories are computed from a two photon symmetrized wave function which has no classical analogue. In 
this sense, the Bohmian trajectories calculated in the following sections represent the first plots of photon trajectories. 

The plan of the paper is as follows. In Section II we study the trajectories in Young's double-slit experiment. In 
section III, we compute the trajectories corresponding to two down-converted photons passing through a double-slit. 
In Section IV we plot the Bohmian trajectories for reflection and refraction through a glass slab. We make some 
concluding remarks in Section V. 



II. SINGLE PHOTON DOUBLE-SLIT INTERFERENCE 



Let us now consider the specific case of double-slit interference of single photons. If the slits A and B have a non- 
zero width d significantly larger than the de Broglie wavelength of the particles (d » A), the slits will convert plane 
incident waves into plane diffracted waves sufficiently far from them (the case of Fraunhoffer diffraction). One can 
see this by carrying out the necessary approximations |f| on the single-particle spherical wave at a point P, arriving 
from a point within a slit at a distance x = ±£ from the origin, and integrating over the slit Q . The wave function 
at a point (x, y) at a sufficient distance D » d 2 /X to the right of the plane of the slits is given by[] 



ip{x,y) = Ma9a 



exp(ifcr J 4) 



Mb9b 



exp(ifcre) 



where qa and gs are the diffraction factors given by 



9a,b 



sin{kyd/2D) 
kyd/2D 

and M.a and Mb are the Kemmer-Duffin wave functions given by 



(20) 



(21) 



M A = 



-E o sin{0 A ) 
-E cos(9 A ) 







B 







(22) 



and 



1 Henceforth we shall write i{j in place of Tip for brevity of notation. 
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/ E sm(d B ) \ 
i -E cos(6 B ) » 





M B = r 





V o / 

where 9 A and 9 B are the angles of diffraction from slits A and B respectively. 

Using the above wave function the components for Bohmian veloccity are given by 

IE B I 

v x = ° ( g A cos(6 A ) + g B cos(0 B ) + g A g B cos[k(r A - r B )](cos(9 A ) + cos(9 B )) 
v y = ( -g A sm(d A ) + g%sin(9 B ) + g A g B cos[k(r A - r B )](sm(0 B ) - sin^)) 

with ip'ip given by 



(23) 



<ty = (El + B 2 Q ){g A + g 2 B ) + 2g A g B ^ 2 cos(^ + B ) + B 2 ^cos[/c(r A - r B )] 



(24) 



(25) 



The Bohmian trajectories for photons can now be plotted for different initial positions along the slits using the 
above expressions for the velocity. We have taken a uniform distribution of the initial positions for both the slits. 
(See Figure 1). The trajectories clearly correspond to the probability density obtained from standard quantum theory 
at any line parallel to the line joining the slits (y-axis). The trajectories are similar to the trajectories of massive 
particles M. 



III. BOHMIAN TRAJECTORIES OF A PAIR OF DOWN-CONVERTED PHOTONS 

Before we proceed to compute the Bohmian trajectories for a pair of photons, the following point needs to be 
clarified. Let us define the rank-2 tensor current 

^(si.ara) = aP(x 1 ,x 2 )^^ ) +$ ) pW -ff^fl^W + ^ a) ^ J - sJa^^.a,) (26) 

for wave functions which satisfy the symmetry i/)(xt,X2) — ij)(x2,Xi). Then the i-th component of the Bohmian 
velocity for the n-th particle (n — 1,2) is 

vl n \ Xl ,x 2 )=c^^\ (27) 

Using similar arguments to those presented in section I, it is clear that this Bohmian velocity is also time-like. The 
expression (27) however appears to be non-covariant because the two sides transform differently. Nevertheless, it is 
possible to write it in a manifestly covariant form by introducing a foliation of spacetime with spacelike hypersurfaces 
£ with future oriented unit normals rj^ix) at every point x of £ such that r]^(x)rj^(x) = 1. Then 

(1)/ \ s ili {xi,X2)rf{x2) 

v\ '{xi,x 2 ) = c — 

(xi , x 2 )?7 M [xi )rf [X2 ) 

V (2) (X1 X2) = C *M(S1.*2W(S1) . m 

1 [U 2) S I1U {X 1 ,X2)^{X 1 )^{X2) (28) 

The fact that EPR entangled states can be written in a manifestly covariant form using this technique of spacetime 
foliation was first shown by Ghose and Home M . The same technique was used by Durr et al. M and Holland M in 
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the context of Bohmian velocities for multiparticle entangled states to demonstrate their relativistic covariance. For 
further details, see P,^p[. 

We now consider an experiment in which a pair of down converted photons is made to pass through two identical 
slits. We will compute the Bohmian trajectories for this case in the limit of Fraunhoffer diffraction. The two-particle 
wave function (in the Fraunhoffer limit, i.e., x\ — x 2 = D >> d 2 /A) is given by 

cxv>(2ik,D) ( i 

i>(VuV2) = ^2 ~ d2 9i92 M A iM B2 exp(-ika( yi - y 2 )/D) + M A2 Mmexp(ika( yi - y 2 )/D) I (29) 

After substituting the expressions for the Kemmer-Dufhn matrix elements and the diffraction factors, we obtain the 
following expressions for the Bohmian velocities of the two photons: 

vi x = {-2 ^cos(^) + glcos(9 B ^j + g\g\ (l + cos(^ + 9 B )) (cos(6 A ) + cos(0 B )) ^ (30) 

V2x = V\ x (31) 

2^sin(0 A ) _ glsm{9 B ^j + g\g\(\ + cos(0 A + B j) (-sin(^) + sin(0 B ))^ (32) 

v 2y = —Vi y (33) 



Vly 



where ip'ip is given by 



{cos(9 A + 9 B ) + 1 



,2 



1 + V ; ' cos(2ka( yi - y 2 )/D) 



(34) 



The second cosine term represents a fourth order interference in the joint detection probability of the two photons 

The Bohmian trajectories are plotted in Figure 2. It can be checked that they agree with the joint detection 
probability amplitude obtained on a plane parallel to the plane of the slits. Again, they are similar to the trajectories 
one obtains for the symmetrized wave function of two massive particles Note that the trajectories are symmetric 
about the x-axis, and the trajectories in the upper and lower half-planes do not cross. 

IV. REFLECTION AND REFRACTION THROUGH A GLASS SLAB 

Finally, let us consider the example of refraction of light through a glass slab. We consider both the air-glass 
interfaces separately and combine the solutions. To obtain the solution of the Kemmer-Dufhn equation in this case, 
we must first solve Maxwell's equations for this case. Let the electric field be polarized along the y direction and let 
the wave propagate along the x direction. The air-glass interface is taken at x = 0. Let the amplitude of the electric 
field be represented by a gaussian wave packet^ centered at xq, i.e., 



E x = £; exp | 



— (x — ct — Xq) 



2 



2(70 



E y =0 



E z = (35) 
Taking into account the boundary conditions at the air-glass interface, one obtains 



2 Such wave packets are nowadays routinely produced in the laboratory in down-conversion experiments. See, for example jo 
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E x (x,t) = 



i?oexp 



~(x — Ct — Xq)' 

2^0 



ffi \Eoexp 



-{x + ct — Xq) 



V 



^ exp( -(--f-°) 2 )U>0 



x < 



(36) 



where n is the refractive index. The corresponding magnetic field is given by 



B z (x,t) = 



L exp 



— {x — ct—Xf)Y 
2^0 



1-n En 
1+n I c 



exp 



V 



-(nx—ct—xo)* 
23^ 



— (x-\-ct— xo) 
2^o 



, X > 



X < 



(37) 



The Kemmer-Duffin-Harish Chandra wave function is therefore given by 

-Ey 

-E z 

Bj; 

B u 
B z 







V o / 



1> 



(38) 



Using these expressions for the electric and magnetic fields, one obtains the Bohmian velocity to be 



/ GXp 



"0 



-{x + at + xgy 
"0 
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cxp 



|, X < 



V 



(39) 



Similarly the solutions for the electric and magnetic fields, and the corresponding expressions for Bohmian velocity 
can be obtained for reflection and refraction at the next glass-air interface placed at x = 0.2. These two sets of 
solutions are combined to obtain the trajectories for photons reflected and transmitted through the glass slab. The 
trjectories for a particular set of initial positions are plotted in Figure 3. 



V. CONCLUSIONS 



Bohm and his coworkers have all along emphasized a fundamental difference between fermions and bosons in that 
fermions, in their view, are particles, whereas bosons are fields. This asymmetry in the Bohmian picture of fermions 
and bosons arose due to the absence, in their view, of a consistent relativistic quantum mechanics of bosons with 
a conserved four-vector current which is time-like and whose time component is positive. Such a formulation was 
provided by Ghose et al. [0,0 and it was shown that Bohmian trajectories for relativistic bosons could be defined Jl3| . 
Just as the actual plotting of Bohmian trajectories for nonrelativistic particles was an important advance jl4j, it is 
equally important to demonstrate the actual nature of Bohmian trajectories for relativistic bosons in simple physical 
situations, particularly because such trajectories were thought not to exist by Bohm himself. This does not in any way 
detract from the significance of Bohm's general point of view regarding the causal interpretation. In our view these 
trajectories constitute a significant support of Bohm's causal interpretation by removing an unnecessary asymmetry 
between fermions and bosons from it. In case there is any truth in supersymmetry, such an asymmetry would be fatal 
for Bohmian mechanics. 
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Figure Captions: 

Figure 1. Bohmian trajectories of photons for self- interference through a pair of identical slits centered at y = —0.0002 
and y = 0.0002. 

Figure 2. Bohmian trajectories for a pair of photons passing through two identical slits. Note that there is no crossing 
of trajectories between the upper and lower half planes. 

Figure 3. Bohmian trajectories for photons passing through a glass slab placed at < x < 0.2. Reflection and 
refraction are seen at both the air-glass interfaces. 
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